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Abstract. In this paper we investigate the use of path- integral formalism and the 
concepts of entropy and trafhc in the context of molecular motors. We show that 
together with time-reversal symmetry breaking arguments one can find bounds on 
efficiencies of such motors. To clarify this techinque we use it on one specific model 
to find both the thermodynamic and the Stokes efficiencies, although the arguments 
themselves are more general and can be used on a wide class of models. We also show 
that by considering the molecular motor as a ratchet, one can find additional bounds 
on the thermodynamic efficiency. 



1. Introduction 

Understanding the functioning of living organisms has always been one of the most chal- 
lenging and most interesting aspects of science. It has been known for a while now that 
molecular motors play an important role in the working of these organisms [H [2] . The 
term is actually used for a whole collection of small machines. However, they all have in 
common that they transform chemical energy into work, and that thermal fluctuations 
are essential to their working. The energy needed comes mainly from the hydrolysis 
reaction of ATP ^ ADP + P, see e.g. [21 H [5]. 

The chemical reaction cycle that drives the motor, is itself driven by maintaining a 
high concentration of ATP in the motor's environment, and is thus essentially out 
of equilibrium. Furthermore, the fact that molecular motors are very small, means 
that we are dealing with mesoscopic systems for which the dynamics involve stochastic 
elements. Computations of e.g. efficiency of these molecular motors must therefore be 
quite different from their thermodynamic analogues. The question of efficiency is indeed 
important and has been widely discussed already in the literature, see [61171 ISllQl fTOl [TT] . 
In this paper we present a general method for obtaining bounds of efficiencies based on 
the relation between time-symmetry breaking and entropy production. These results 
are quite general and do not depend specifically on a model, even though we will il- 
lustrate the argument by using one specific model. We give bounds on both the usual 
thermodynamic and the Stokes efficiency. Molecular motors are part of a wider class 
of systems called Brownian motors or ratchets [121 [131 (H] We will use this relation be- 
tween molecular motors and ratchet dynamics to derive other, less explicit, but possibly 
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better bounds for the efficiency. Also here we show that general probabilistic arguments 
can replace the thermodynamic procedure of computing efficiencies once there is a clear 
identification of the path-dependent entropy fluxes. 

A well-known example of a molecular motor is kinesin, which moves cargo inside the cell 
away from the nucleus by "walking" along microtubules, with a step size of 8nm and 
loads of the order of pN, multiplying to work of the order of fc^T at room temperature 
[3]. This is a small scale for which a stochastic description is best fltting. 
The hydrolysis reaction mentioned above will be modelled by a chemical cycle, described 
by a Markov jump process. In a molecular motor this chemical cycle induces movement 
of a load in a spatial direction, thus doing work. The spatial movement of the system 
will be described here by an overdamped Langevin equation, due to the viscosity of the 
medium in which these motors live. The coupling between the chemical cycle and the 
spatial movement is a so-called "ratchet effect:" to each chemical state of the system 
corresponds a different potential landscape, and by switching between these, a directed 
movement arises. The system is driven out of equilibrium. In the chemical state this 
is due to the maintained chemical potential, which is the actual energy input. The 
spatial component is in a sense also driven out of equilibrium by a non-gradient "load" . 
However, the system will move in a direction opposite to this load, performing work as 
an output. 

We start in section [2] by describing the model and introducing some preliminary con- 
cepts. In section |3] we will look at the entropy flux, the entropy production rate and 
the positivity of the entropy production using the path space formalism. In section H] 
it will be shown that, in accordance with the second law of thermodynamics, both the 
thermodynamic efficiency and the Stokes efficiency are bounded. Section [5] treats the 
motor as a ratchet. Using symmetry arguments from the ratchet formalism, we examine 
the efficiency more closely. 

This paper is intended to illustrate two things: first of all a method to check the usual 
bounds on efficiencies - that they are smaller than 1. The model and the bounds on 
its efficiencies have been described already in |6] . We focus here on the physical (model 
independent) argument, which revolves around the path-space formalism, entropy pro- 
duction and the second law of thermodynamics. Secondly, we illustrate how the ratchet 
formalism can produce possibly better bounds. In the model we use to illustrate the 
method, the new bound is given by fl29|) . 



2. The model 

We consider a model for a molecular motor, coupled to a thermal bath at inverse temper- 
ature /3, which has 2 degrees of freedom : a position x on a one-dimensional circle with 
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length L, modelling a periodic potential landscape, and a chemical state i G {1, . . . , A^} 
with periodic boundary conditions: = 1. As time evolves, the system will therefore 

cover a trajectory (path) in the space of possible configurations. Such a path starting 
at time zero and ending at a later time T will be denoted by c^; = {xt,it)o<t<T- 



dxt 



dt 



'2DdB, 



(1) 



The dynamics in the spatial direction x is taken to be an overdamped diffusion 
process (defined with the Ito convention): 

where ^ is the viscous drag coefficient, D is the free diffusion coefficient, related by the 
Einstein relation to the inverse temperature: (3 = jj^. Furthermore / is an external 
load, which is a nonconservative forcing, and 0o(2;t) ^t) is a potential, which depends on 
the chemical state the motor is in. The prime denotes here the derivative with respect 
to X. Finally, dBt is standard Gaussian white noise. 



The dynamics of the chemical state is governed by the transition rates kx{i,i + 1) and 
kx{i,i — 1), which depend on the position on the circle. Here we impose the physical 
condition of local detailed balance: 

kx{i,l + 1) ^ g/3[0o(x,i)-0o(a;,i+l)]+/3A/x(i,i+l) 

kxii + l,i) ^ ' 

which just states that the rate for jumping to the right divided by the rate for jumping 
to the left is equal to the exponential of the entropy produced by jumping to the right. 
Indeed, the entropy change in going from state z to z + 1 at position xis (3 times the heat 
flux between system and environment, which itself is equal to minus the change in energy 
(0o) plus the chemical work A/x(i, i+1). Note that A/i(2, i+l) = — A/i(2+l, i). In general 
every full chemical cycle (z — > i + A^) consumes a chemical energy A/i = AyU,(i, i + 
When this A/x is not zero, it is seen in ([2]) that detailed balance is broken. The chemical 
potential thus provides the energy input for the motor, driving the chemical reaction 
cycle. 

Let pt(x, i) be the probability density of finding the motor in position x in chemical 
state i at time t. From the dynamics described above, one can prove that this probability 
density evolves according to a set of coupled Fokker-Planck equations : 

dpt{x,i) _ dJp,{x,i) 



. ^Pt(a;,« - 1,0 - ^Pt(a;,^,^ + 1) (3) 

Where JpXx.,i) is the probability current in the spatial direction for the i-th state and 
Jpj(a;, i,i + 1) is the probability current from state i to state i + 
These probability currents are given by 

/pt(x, z, z + 1) = pt{x, i)kx{i, i + 1) - Pt{x, i + l)kx{i + 1, i) (4) 
JpXx, i) = ^Ptix, i) [f{x) - 0o(x, i)] - Dp[{x, i) (5) 
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Taking the left-hand-side of Q zero gives the equation for the stationary measure, which 
we will henceforth denote by p(a:, i). It is easily seen that when / and the A/i(z, « + 1) are 
put to zero, the process satisfies detailed balance, with the equilibrium measure given 
by the Gibbsian distribution 

In the following it will be convenient to define such an equilibrium process as a reference 
process. More specifically with rates kP^^ given by 

l3Afi(i.i+l) 



„ /3A)j(i,i+l) 

kl{i + 1,0 = k^{i + 1, i)e ^ 
and with the same Langevin dynamics as in ([T]) with / = 0. 

3. Entropy Production 

3. 1 . Path space formalism 

A useful way of describing a nonequilibrium system is via the distribution over 
its possible trajectories (paths), denoted by Pp^iuj) with po some initial measure. 
A consequence of the Markov property is that one can write that Pp(,(u;) = 
PQ{xQ,iQ)P{ijj\xQ,iQ), where the second factor is now independent of the initial measure. 
Such path-space distributions are best described as relative densities with respect to a 
reference process, denoted by P°. For this purpose we take the reference equilibrium 
process as defined at the end of the last section. Given an initial distribution p^^x^i), 
the ratio of the probality of a path uo is given by 

dPpoiuj) ^ po(xo,zo) P{Lu\xo,io) ^_ Po(xo, ip) 

dP^^{uj) pQ{xo,iQ) P^{u\xo,io) ' pQ{xo,io) 

where A is called the 'action.' The action will be a sum of the action for a jump process 
(see e.g. Appendix 1 of [15]) and the action for a diffusion process (see e.g. [T6]): 

A{uj)= f dxtof{xt) + ^^ f dtf\xt) 
^ Jo Jo 

-|f*/K)[/w-2^;K>.)i-Eiog(M) 

+ r dtY,[k.Ah,j)~kl{H,j)] (8) 
Jo ^ 

where the first term of the action is a Stratonovitch-type stochastic integral, and Y2t<T 
means the sum over all the times that the system jumps to another chemical state. By 
it- and it we denote respectively the chemical state before the jump and after the jump. 
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With the path space distribution, one can compute expectation values of 
observables. When started from an initial measure Pq, the expectation value of an 
observable O^u) is denoted by 

{0{u)l^ := / dPMO{u) 



where the integral means the integral over all possible paths. 
3.2. Entropy flux 

It is well known ([17]) that the time- antisymmetric part of the action of a stochastic 
process gives the entropy flux S{u) between the system and its surroundings. Let 9 be 
the time reversal operator, acting on paths as = {xx-t, iT-t)t- The entropy flux is 
then given by 

S{uj) = A{euj) - A{uj) 



(3 r dxtof+Y. 

J +^1^ 



t<T 

fT 

t<T 



I dxtof + pJ2^^^i't-,^t) (9) 

Jo 



The two terms in this path-dependent entropy production already have a clear meaning. 
The first term comes from the motion in the spatial direction. In fact this first term 
is the entropy flux one would see, if the motor was only allowed to move in the spatial 
direction (for example by setting all rates k to zero). Vice versa, the same is true for 
the second term, which comes from the motion in the chemical direction. We see that 
one full chemical cycle gives an entropy exchange between system and environment of 
PAfi, as expected. 

The path space average of the entropy flux, starting with an initial distribution p, can 
be computed, again through the theory of Markov jump processes and diffusions: 

('5(^))po = /9 / dt [ c/x V Jp,(x,i)/(x) 
Jo Jo j 

+ p[ dt I dxy'/p,(x,z,z + l)A;u(i,z + 1) (10) 
Jo Jo j 

3.3. Entropy production rate 

The total entropy production of the universe Su caused by this process is not only 
the entropy flux, but also the entropy change of the system itself. The entropy of the 
system is given by the relative entropy of the measure it is in, relative with respect to 
the equilibrium reference measure: 

/ dxp{x,i)\og ' (11) 
,_i ^0 Po[x,t) 



Entropy and efficiency of a molecular motor model 



6 



The change in the entropy of the system is then given by 

siPr) - sipo) (12) 

Summing the two contributions ( (fTOi) and (fT2l) ) to the total entropy change Su, it is 
then most elegantly written as 

Sl(P.)--(lo,'-^) (13) 

where pt is the time-evolved measure determined by the Fokker-Planck equations ([3]). 
The Markov property of our process makes that we can write the entropy production 
as 5*^(^0) = /o dt a{pt) where cr(pt) = lim 7^S'J~*(pi) is called the entropy production 

rate, straightforwardly computed to be 

<V) = K I dxY,^-^ (14) 
J ^ P[x,i) 

+ (3 / dx^ Ip{x, + [A/i(2, i + + ^{x, i) - <l>(x, i + 1)] 

i 

where $p(x, i) = (t>o{x, i) + log p{x, i). 

Just as for the entropy flux, this entropy production rate consists of two separate terms, 
each having a clear meaning. For the first term, think of a new dynamics in which 
all chemical transition rates are set to zero. Then there is only motion in the spatial 
direction, and the entropy production for this process would exactly be the first term 
in (fT4l) . Suppose instead that we forbid all motion in the spatial direction, by letting 
^ ^ 00. Then only the chemical reactions determine the process, and the entropy 
production for such a process is exactly the second term in f|T^ . We will therefore write 

<p) = (yD{p) + (y.j{p) (15) 

to denote that the entropy production is the sum of the entropy production of the 
diffusion process and of the jump process. 

3.4- Positivity of the entropy production 

The definitions of entopy flux and entropy production rate through path space measures 
have the advantage that one can easily prove their positivity: because of normalization 
of the path space measure we have for any measures q and time T, and actually for 
any (well-behaved) dynamics: 

MSSA)\ =(^M^\ =1 (16) 

Jensen's inequality then dictates that 

U^X) >0 (17) 



Using the definition of the entropy production rate a we find 
a{p) > 
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Which is in essence the second law of thermodynamics: change of total entropy of the 
world is positive. As said before, this entropy production rate can be written as a sum 
of an entropy production rate due to the diffusion in the spatial direction (ctd) and an 
entropy production rate due to the jump process in the chemical direction (cj). As 
the argument of positivity above is quite general and does not depend on the specific 
dynamics of the system, it is also valid for dynamics where we set some rates to zero or 
let ^ ^ oo. In particular this means that the two terms of ([T5l) are separately positive 

^^(p) > 



4. Efficiency 

To prove that the thermodynamic efficiency is bounded by one is not difficult. One just 
has to realize that the total stationary entropy production of the process is the heat flux 
between the system and the environment. This heat flux consists of the energy going 
into the system Qm, in our case through the chemical reaction cycle, and the power 
output Wout by moving a load. The fact that the stationary entropy production has 
to be positive then immediately gives that Qin — Wout > 0, or equivalently that the 
efficiency ^yf^ < 1. 

But there is more. First of all, one has to realize that molecular motors like the one 
described here, are essentially two processes which are coupled. The entropy production 
of the whole process is then the sum of the entropy productions of the two separate 
processes. This is just a consequence of the Markov property. There is however one 
important difference: the stationary measure of the whole process is not the product 
of the stationary measures of the separate processes. This has as a consequence that 
the stationary entropy production of the whole process is rather the sum of transient 
entropy productions of the separate processes. Still, these entropy productions are both 
positive, which leads to the bound on the Stokes efficiency, as illustrated for our model 
below. 



4.1. Thermodynamic efficiency 

Here we will work in the stationary case, i.e. the measure will be chosen to he pt = p 
the stationary measure. There is an immediate consequence of the (stationary) Fokker- 
Planck equation ([3]): 
dV 

— =0, with y:=L^J,(x,2) 

i 

z/(z) — z/(z + 1) = 0, with u := dx Ip{x,i,i + 1) 



which unambiguously defines the average velocity V and the mean chemical reaction 
rate i'. As usual, the efficiency is the (average) mechanical work that the motor can 
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perform on the load, divided by the (average) chemical energy that is provided. With 
the above definitions, the thermodynamic efficiency is then clearly equal to 

where f = j; Jq dxf{x) is the average force on the circle. The minus sign in the definition 
is there because for the motor to work, one needs that uAfi > and fV < 0. As we 
know, the second law of thermodynamics should give a bound on this efficiency, namely 
that it is smaller than 1. Indeed, in the stationary case the entropy fiux ( fTOl) becomes 

{S{uj))^ = (3T[Vf + uAf,] (20) 

When the system is in the stationary state, its Shannon entropy (fTTl) is constant in time, 
so that 

{S{u))^ = Ta{p)>0 (21) 
Together with ( l20l) . this gives that rjT < I. 

4-2. Stokes efficiency 

In many cases the goal of a molecular motor is to carry cargo across a viscous 
environment (e.g. protein transport within a cell). In this case the useful power of 
the motor is just the viscous drag force times the average velocity: C,V ■ V. While the 
total power input from the environment is the chemical input through the chemical 
potential, and mechanical input through the external force /. Note that / is now not 
regarded as a load, but as an external force. Because of this, the Stokes efficiency 
was introduced in ([IB]) to have a more useful definition of efficiency than the previous 
thermodynamic one: 

We will show here that for this model, the Stokes efficiency is bounded by one, just 
as the thermodynamic efficiency. Again, this will follow from the second law of 
thermodynamics, by considering the entropy production rate o"(/i) as the sum of the 
two entropy production rates ao and cxj for the mechanical and the chemical cycle. We 
already proved that these two quantities are positive (see section [3^ . Note however, 
that these two quantities are entropy productions for restricted dynamics, one where 
the chemical movement is prohibited and the other where mechanical movement is not 
possible. This also means that the stationary measure of the molecular motor is not 
the stationary measure for these restricted processes. Luckily, this poses no problem, 
because the positivity of the entropy was proven for any measure p. We thus have that 
CToip) > and aj{p) > 0. We find then 

MP) ^^J^-?'"-')(^)^ 
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This inequality leads to the bound on the Stokes efficiency: 



Vs 



N 



f Jp{x)dx 

i=l 



fextV + uAjJ fextV + uAjJ 



< 



(3 [UtV + i/A/i] 
< 1 

where the ffist inequality is an immediate consequence of the positivity of the variance 

of ^. 

p 



5. The motor as a ratchet 

In this section we will show that the molecular motor as modelled here is a exactly a 
ratchet as defined in [12]. We will then use properties of ratchets to derive another, 
possibly better, bound on the thermodynamic efficiency 



5.1. Another symmetry 

The quantity that determines the power output of the motor, is the current in the spatial 
direction. For a specific trajectory u, this current Jt is defined by the following relation: 

j dxg{x)jT{x) := dxtog{xt) (24) 

where again a stochastic integral is taken with the Stratonovich interpretation. With 
this definition, one can compute that {jrix))^ = j"^ dtJ^^{x). 

It is obvious that jt is independent of the specific sequence of chemical states it went 
through, it only depends on the movements in the spatial direction. It will therefore be 
convenient to define a new symmetry operation F as follows: Tu := {xt,N + 1 — it)t- 
The current jrix) then satisfies 

^Jt{x) = Jt{x) 
dhix) = -jrix) 

which makes it a ratchet current, exactly as defined in [T2] . 



We now have two symmetry operations: the time-reversal 6, and the "chemical 
reversal" F. With this, we can repeat the argument of the positivity of the entropy 
production, but now with our new symmetry operator. Let us therefore split the action 
([8]) into four parts, according to these symmetries. The antisymmetric part of the action 
under time reversal was already seen to be the entropy production. The time-symmetric 
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part of the action was also earlier examined in e.g. [T^l [ISj, and named traffic T, so 
that 

A = l[T-S] 

Using the symmetry operator F, the action splits further as follows: 

where the plus and minus denote respectively the symmetric and antisymmetric parts 
of the functions under F, e.g. = [TS — S]. For the argument to work, it is necessary 
that the reference equilibrium process is also symmetric under F. Even though it is 
always possible to find such a reference, to make computations simpler, we will assume 
that the reference process we defined earlier is already symmetric under F: 

4>o{x,i) = 4>o{x,N + 1 -i) 
kl{i,j) =A;°(iV + l-z,iV + l-j) 

With this, we repeat the argument of section 13. 4[ but now for TO instead of 6. This 
argument gives here that 

1 p{xo,io) \ 1/c / ^ , u /i dPf,{uj) \ . „ 

p{xT,N + 1 -zt) / p 2 P \ dPpiVeuj)/^ 

which means in particular that 

As this equality is true for any time T, and as the ffist term on the right-hand side of 
this inequality is not extensive in time, one sees that 

^(T_H)^>-2/3V/ (26) 

We see therefore that the time-symmetric part of the action plays an important role 
in ratchets and thus in molecular motors. Let us therefore try to compute 71 for our 
model. 

5.2. The efficiency revisited 



In this case it is easily seen that 
TJuj 



2/ dty^[KXN + i-tu3)-KAh:3)] 

Jo ^ 

After a straightforward calculation, one finds 

{T4u))^ = 2TJ2 cosh('^A/i_(^,^ + l)') -1 r(^,z + l) 

+ 2uT sinh ^^A/i_(2, i + 1)^ 



(27) 
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where we have defined A/i_(i, i + 1) = Afi{N + 1 — i,N 



A/i(i, i + 1), and 



r(z,z + 11 



dx p{x, i)kx{i, i + 1) + p{x, i + l)kx{i + 1, i) 



is the symmetric counterpart of the current in the chemical direction. 
It follows now from (l26l) that 



cosh ( ^A/i_ 



(z,z + 1) 



-f3Vf< 

i 

+ V ^ sinh A/i_ (i, i + 



1 



r(z,z + 1) 



1) 



We find a new boundary for the thermodynamic efficiency 



< 



T{i,i + 1] 



cosh ( ^Afi_{i, i + I] 



+ 



-J— sinh f 



/5 



A/i_(z,z + 1) 



(2^ 



(29) 



Even though this is not a very explicit formula, one can see the importance of 
fiuctuations: the quantity r(i,i + 1) deals with the average number of jumps between 
states i and i + 1, thus being a measure of the total activity in the chemical direction. 
When r(z, i + l)/^ is big, then the efficiency of the motor may also be big. Furthermore, 
this traffic plays a fundamental role in fluctuation theory, see e.g. [ISl [IS]- 



6. Conclusions 

We showed in this paper that the boundedness of the efficiencies (both thermodynamic 
and Stokes) for a large class of models for molecular motors can be deduced with model- 
independent arguments. These general probabilistic arguments are based on the path- 
space formalism and the relation between time-symmetry breaking and entropy produc- 
tion. By also showing that such a motor can be seen as a ratchet other boundaries were 
found for the thermodynamic efficiency. 



To find more explicit results approximations should be made, or otherwise a 
numerical analysis. This is beyond the scope of this text, e.g. in [HI [10], [11] such 
approximations and numerical analyses are made, however by using a different method. 
Therefore it could be interesting for further research to make these approximations 
and/or numerical analyses, starting from our method, to compare with the results in 
these papers. 
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